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МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ШУМООБРАЗОВАНИЯ ВИБРО- 
УДАРНОГО УПРОЧНЕНИЯ ЛОНЖЕРОНОВ ВЕРТОЛЕТОВ 


Приведена математическая модель шумообразования стенда виброударного упрочне- 
ния лонжеронов вертолетов. Энергетическими методами рассчитываются скорости 
колебаний отдельных подсистем несущей системы в целом. Получены зависимости 
вводимой от рабочей среды вибромощности в стенки контейнеров, что существенно 
уточняет теорию расчета шума подобного оборудования. 

Ключевые слова: шумообразование, виброударное упрочнение, лонжероны верто- 
летов. 


Введение. Упрочнение внутренних и наружных поверхностей лонжеронов в 
условиях ОАО «Роствертол» производится на однокоординатных станках. 
При всех преимуществах этого оборудования с точки зрения технологии 
упрочнения оно имеет очень серьезный недостаток: создает чрезвычайно 
высокие уровни шума, на 30-33 дБ превышающие предельно допустимые 
значения. Для обоснования выбора варианта шумозащиты при проектирова- 
нии подобного оборудования необходимо расчетным путем определить ожи- 
даемые нормы шума. 
Исследование. С акустической точки зрения каждый элемент вибростенда 
является конструкцией энергетически замкнутой, с небольшим коэффициен- 
том потерь колебательной энергии. В качестве моделей источников шума 
для элементов вибростенда принята ограниченная пластина, для которой 
звуковая мощность определяется по формуле [1]: 
рось5 7 
р! (1) 
10 
где ро и а - плотность воздуха (кг/м?) и скорость звука в воздухе (м/с); 
$ — площадь пластины, м7; И; — скорости колебаний пластины на соб- 
ственных модах, м/с. 
Звуковая мощность оборудования как системы одновременно звуча- 
щих источников определяется по принципу энергетического суммирования и 
имеет вид: 


[= 1915 


т 0,1 
Е. =10 10° +107, 2) 
1 
где В — звуковая мощность, излучаемая контейнером, дБ; Г. - звуко- 


вая мощность, излучаемая рамой, дБ. 

Определение звуковой мощности вибростенда сводится к расчету 
скоростей колебаний рамы и контейнеров на собственных формах колеба- 
ний. В такой постановке расчет вибраций элементов вибростенда на основе 
теории колебаний пластин — задача чрезвычайно сложная. Поэтому для 
практических расчетов виброскоростей элементов стенда использованы 
энергетические методы расчета, которые традиционно используются для 
акустических расчетов энергетически замкнутых систем. Их применение наи- 


217 


Технические науки 








более оправдано для базовых деталей коробчатой формы, что и соответ- 
ствует рассматриваемому оборудованию: 
6 6 


и + 0 Е ю 0 ны г 0,5, ы 0 ' (3) 
1 1 
где И/- вводимая в раму вибромощность, Вт; а, и („; - коэффициенты 
передачи вибромощности между стенками рамы; 4“„ - коэффициент 
поглощения в соответствующем элементе вибростенда, 1/м; 4;„- по- 
токи вибрационной мощности, Вт/м. 
Для стальной конструкции коэффициенты поглощения и потоки 
вибрационной мощности (согласно работе [2]) приведены к виду: 
0,5 
Е в 
$. = 0,0317 =, (4) 
й, 
где 1 и — частотно-зависимые коэффициенты потерь колебательной энер- 
гии; К — собственные частоты колебаний стенки, Гц; №, — толщина 
стенки, м; 


Ч: = 200 И Йа (5) 


где 77%  - распределенная масса стенки, кг/м?. 


Коэффициенты передачи вибрационной мощности между элемента- 
ми рамы определяются следующим образом [2]: 
от. 
к = Ши] и (6) 
Л 
Для Г-образного соединения пластин, что и соответствует конструк- 
ции вибростенда 


-2 
= ий й, при а» = Й, 


р р. 
| Г = Г (7) 
11 — \ 
при < У} = 5 (1+н,) Е 5 | Е 5 | 2 
11 м т — 
при ы 1 Уи, = Е (1+ н,) Е 5 ы + [АЕ 5 ы 2 


Еи Е- полные эллиптические интегралы, таблицы которых приведены в 
работе [2]. 

Поскольку возбуждение вибраций рамы создается вибраторами, то 
вводимая в основание рамы вибромощность задается согласно выражению 


Й’= 0,0023Ма^и*, (8) 


где М -— приведенная масса вибратора, кг; а — амплитуда колебания, м; 
п — частота вращения вала вибратора, об/мин. 
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С учетом этих зависимостей система уравнений энергетического ба- 
ланса рамы имеет вид: 
2,5 
р ВИ | 
-5 23 Г 0,5 2 
5+10`_ Мат ны =" 1, (Л) ту, 
2,5 2,5 р 7 
в [725+ 12°) № 


у 5 (11]” р | 
р 18,81 „Лит, Р. . . 25\2 УЕ й,/, т, 7, $09) 
г [127+ 12°) №, 


Из системы определяются скорости колебаний каждой стенки, а на 
их основе и уровни звуковой мощности. Аналогичным образом определяют- 
ся и скорости колебаний стенок контейнеров. Отличие их расчета от рамы 
заключается в задании вибрационной мощности. Мощность, вводимая в 
элементы корпуса контейнера, складывается из мощности, передаваемой 
от рамы (с учетом ослабления на ложементах), и мощности от воздействия 
рабочей среды. Кроме этого учитывается воздействие рабочей среды на 
боковые стенки и основание контейнеров. Формулу мощности, вводимой в 
корпус контейнеров, представим в виде [99]: 





у 
С в 

Й’. — и  ., (10) 
о 

где Г и Г - виброскорости колебаний нижней и верхней пластин 


[9 в 
рамы, м/с; И/„: — вводимая в контейнер вибромощность от рабочей 
среды, Вт. 

Внутри вибрационного контейнера выделена некоторая поверх- 
ность () , подвергающаяся воздействию частиц технологической среды. 
Это может быть как поверхность любой стенки контейнера, так и поверх- 
ность упрочняемой детали (для дальнейшего рассмотрения это не имеет 
значения). Площадь поверхности ©) обозначим $. Основной характери- 
стикой виброакустической активности действия среды на данную поверх- 
ность является: 


Й’ = (Е и’ у 2 (11) 


где Е -— сила воздействия; и” — относительная скорость падения среды на 
движущуюся поверхность, индекс ©) при угловых скобках означает сум- 
мирование по всей поверхности () ‚ а & - достаточно малый промежуток 
времени, за который производится усреднение. Величина вводимой в по- 
верхность () мощности И/ является функцией времени, которая, учитывая 
периодический характер вибрационного возбуждения среды и физический 
смысл мощности, должна быть также периодической и положительной. 

Для определения величины И/ в соответствии с (11) рассмотрим 
схему, представленную на рис.1. Выделим на поверхности 0 бесконечно 
малый элемент площадью сб. Направим из центра площадки а внутрь сре- 
ды единичный вектор нормали п, который вместе с парным ему единичным 
касательным вектором Т образует локальный ортогональный репер, дви- 
жущийся по заданному закону вместе с площадкой &$ со скоростью ЦБ. 
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Пользуясь  соотноше- 
нием между тензором напря- 
жений О в среде и силой 
воздействия этой среды на 
ограничивающую поверхность 

Е= -с. п: (5, 

(12) 





а также связью между скоро- 
стью Ч среды в неподвижной 
системе и скоростью ча” в ло- 
кальной системе, движущейся 


вместе с репером (п, т) 
Рис. 1. Расчетная схема для определения интегральной мощности воздействия гра- 


и =9-О,, нулированной среды на поверхность 


получим 
и: (Е и’) =- (6 п) (-0,)) 4). (3) 


Величина, стоящая под знаком усреднения и являющаяся результа- 
том скалярного умножения векторов силы и скорости, изменяется вдоль 
всей поверхности ©), причем тензор 0’ и вектор и определяются динами- 
кой среды, векторы п и ЦЬ - геометрией и кинематикой границы. 

Рассмотрим подробнее величину б. и способ ее определения по 
фазовой траектории гранулярной системы. Произведение суммы таких из- 
менений импульса для всех частиц, столкнувшихся с выбранным сплайном 
за время моделирования &Ё„ на вектор п, отнесенное ко времени &, и ве- 
личине площадки 45 , даст искомую величину тензора напряжений: 


==“ у [лип ). (14) 


11° 45 ми) 





В соотношении (14) величина №. (45) обозначает число столкно- 
вений частиц среды со сплайном 45. Пусть площадка (5 имеет длину 
(вдоль которой расположены частицы одного слоя) В и ширину а (вдоль ко- 
торой размещены слои, совершающие плоское движение). 

Тогда ширина площадки может быть выражена через число слоев 
следующим соотношением: 

нЕ 3 
ее (15) 
2 К» 4 

Чтобы учесть вклад в величину напряжений движения К слоев 
необходимо умножить соотношение (21) на число слоев А. Коэффициент, 
стоящий перед знаком суммирования, будет равен: 

тк тк _ т.К _ 14 т 16 
„45 п аь Е.Б. К. 3/4 \3 в -Ь.6 ° о 

Отметим, что величина, даваемая соотношением (16), определяет 
тензор напряжений в любой системе координат, в том числе и движущей- 
ся, порождаемой площадкой (5. Учтем при этом, что изменение импульса 





220 


Вестник ДГТУ, 2009. Т9. №2(41) 








каждой столкнувшейся С границей частицы также инвВариантно относитель- 
НО выбора системы отсчета: 


Ди=у-и=(у-0,} -(У-0,} у-. (17) 
Тогда выражение 


4 т 
= ое Д . 


вычисляемое в этой новой системе координат, выделяет нормальную со- 
ставляющую напряжений. Дуальное ему выражение 


14 т 
= Е Д й 
5: 3 1, ы Ь. [6] 2 ы 9 (19) 


дает величину касательных напряжений. Это позволяет преобразовать вы- 
ражение (17) следующим образом: 


4 тэ 
3 г, бЬ №, (4) 


(у- м) (и-(,]}, (20) 
0 

причем индекс © имеет смысл суммирования по всем прямолинейным 

сплайнам модельной поверхности. 

Оценим удельную мощность ($ = 1 м”), вводимую в граничную по- 
верхность при синфазных нормальных ударах с частотой 15 Гц одного 
слоя, размещенных без зазора стальных шариков диаметром 5 мм. Для это- 
го предположим, что длина сплайна равна половине диаметра шарика 
Ь= 4/2, следовательно, за один период сплайн границы в среднем за пе- 


риод получает 1/2 удара. Положим также, что граница покоится (Ч, = 0), 


коэффициент упругого восстановления равен А, удары наносятся с часто- 
той # и за один период слой шаров наносит по поверхности только один 
удар. В условиях этих предположений получим выражение удельной мощ- 
ности: 


И’ = 0,6р „АГ(1+ Ки’, (21) 
где р. - плотность среды, кг/м”; #-— частота вибрации, Гц. 


В соответствии с технологической инструкцией упрочняющая обра- 
ботка изделия производится в 4 перехода по схеме, приведенной на рис.2. 

















Рис.2. Схема обработки лонжерона (в порядке следования 1, 2, 3, 4): 
буквами обозначены поверхности, подвергающиеся воздействию рабочей среды 
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Закон движения контейнера с деталью отвечал режимам работы 
реального технологического оборудования: амплитуда вертикальных коле- 
баний 6 мм, частота 15 Гц. Имитационное моделирование программным 
комплексом СбгапМо$ производилось с шагом по времени 5-10 мкс до дости- 
жения реального времени — одна секунда, т.е. полных 15 периодов 
колебаний. 

Вычисляемым энергетическим параметром была удельная мощ- 
ность, вводимая средой в границу, т.е. мощность, приложенная к 1 кв.м. 
площади поверхности стенки. Для определения средней удельной мощно- 
сти на поверхности стенки ее мгновенное распределение численно инте- 
грировали по всей длине куска стенки: 


1 
М№М()= — о(г,Оаг. (22) 
АР 

Распределения удельных мощностей были построены для поверхно- 
стей, подвергающихся воздействию среды стальных шариков. По этим рас- 
пределениям были вычислены средние за период удельные мощности, а 
умножение последних на Площади соответствующих поверхностей дало 
полные мощности, характеризующие акустическую активность обрабатыва- 

ющей среды. Результаты приведены в таблице. 


Средние мощности воздействия технологической рабочей среды 
на стенки контейнера и обрабатываемого лонжерона 


Мощность воздействия среды 
Площадь, средняя средняя 
м? удельная, полная, 
/ 


№2 
№3 


№4 
1,06 1,09 


Выводы. Зависимости мгновенной удельной мощности от времени показы- 
вают, что взаимодействие среды со стенками детали и контейнера НОСИТ 
импульсный характер, продолжаясь около половины периода колебаний. 
Такой характер изменения мощности должен приводить к генерации широ- 
кополосного спектра шума. 

Приведенные зависимости позволяют определить уровни шума дан- 
ного типа оборудования и учитывают компоновку стенда, технологические 
режимы работы. Для выполнения расчетов ожидаемых уровней шума при 
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проектировании подобных технологических процессов и оборудования, на 
котором они будут выполняться, необходимо определить вводимую в кор- 
пуса контейнеров вибромощность со стороны рабочей среды. 
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